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Abstract 

Polymers adopt extended coil and compact globule states according to the balance between en¬ 
tropy and interaction energies. The transition of a polymer between an extended coil state and 
compact globule state can be induced by changing thermodynamic force such as temperature to al¬ 
ter the energy/entropy balance. Previously, this transition was theoretically studied by taking into 
account the excluded-volume interaction between monomers of a polymer chain using the partition 
function. For binary mixtures of a long polymer and short polymers, the coil-globule transition 
can be induced by changing the concentration of the shorter polymers. Here we investigate the 
transition caused by short polymers by generalizing the partition function of the long polymer to 
include the excluded-volume effect of short polymers. The coil-globule transition is studied as a 
function of the concentration of mixed polymers by systematically varying Flory’s y-parameters. 
We show that the transition is caused by the interplay between the excluded-volume interaction 
and the dispersion state of short polymers in the solvent. We also reveal that the same results can 
be obtained by combining the mixing entropy and elastic energy if the volume of a long polymer 
is properly defined. 

PACS numbers: 


*Electronic address: k-seki@aist.go.jp 


1 



I. INTRODUCTION 


A single polymer molecule can transition between elongated coil and compact globule 
states upon changing the temperature, pH of salts or addition of other polymer molecules. 
In general, the extended (coil) state of a long polymer is stabilized when the interaction 
oetween the polymer segment and solvent is favorable, typically occurring in good solvents. 
^ In contrast, the contracted (globule) state of a long polymer is stabilized when the 
interaction between the polymer segments is favorable, which occurs in poor solvents. By 
varying the solvent quality, the coil-globule transition of a polymer can be induced. 

Transitions between coil and dobule states can also be induced by adding solutes that can 
be regarded as another solvent. |2|-l6| When a good solvent is added to another good solvent, 
a coil-to-globule transition may be observed as if the long polymer is in a homogeneous poor 
solvent. Similarly, an unintuitive transition from globule to coil can be obtained by mixing a 
polymer with poor solvents. The added solute can be short polymers. When short polymers 
are added, both coil-to-globule and globule-to-coil transitions have been observed. 
Inspired by the complex transitions of polymers induced by addition of polymer solvents, 
here we study theoretically the structural transformations of a long polymer induced by 
polymer solvents. From another point of view, addition of polymer solvents puts a polymer 
in crowded environment (macromolecular crowding). Coil-globule transition can be regarded 
as the structural transition of a long polymer caused by “macromolecular crowding”. 

A qualitative description of the coil-globule transition was initially formulated using 


mean held theory. 10l-ll2| Later, the nature of this phase transition, that is, whether it 


was continuous or discontinuous one, was extensively investigated using simulations, ^ 
consistent held theory and sophisticated models that considered polymer stihness 


Q 
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However, the additional ehect of short polymers on the coil-globule transition of long 


polymers has not yet been fully studied. Part of the observed concentration dependences has 
been reproduced by phenomenolo^al mean held theories, but the results were not further 
justihed using statistical theories. [7H9| Although the coil-globule transition of polymers has 
been widely studied both theoretically and experimentally, it remains difficult to construct 
theories when a single molecular transition is induced by a change in the concentration of 
other polymers. 

Here we study the transition between coil and globule states of long polymers induced by 
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addition of short polymers by a statistical theory using the partition function. The excluded- 
volume effect of short polymers is taken into account as a factor that decreases the emptiness 
in the partition function in a mean held sense. We also include Flory’s y-parameters to take 
into account the interactions between the long polymer, short polymers and solvents. From 
the extremum of the partition function, we obtain an optimal expansion factor for a long 
polymer. 

Among many possible conformations, the optimal ones can be determined from both the 
entropy in the long polymer and its mixing entropy with other polymer molecules as well 
as their interaction energy. There is some ambiguity in combining the mixing entropy of 
polymers and elastic energy originating from entropy in a single polymer molecule. Using 
the partition function, such ambiguity is removed. We show that the same equation for 
the expansion factor can be derived by combining the mixing entropy of polymers and the 
elastic energy originating from the entropy in a single polymer molecule. 

In our approach, the optimal expansion factors are obtained from the extrema of the 
partition function, whereas the steady-state expansion factor was obtained by imposing the 
osmotic pressure balance at the boundary of the region occupied by a long polymer in pre¬ 
vious approaches. [t-Q] In these theories, the region occupied by a long polymer is tacitly 
assumed to be in a phase different from that of the outer region for short polymers and 
the solvent. Furthermore, in principle, the pressure balance should include the hydrostatic 
pressure. In these theories, osmotic pressure was taken into account, but the hydrostatic 
pressure was not fully described. As in the previous theories, we ignore the hydrostatic pres¬ 
sure. Our results, nevertheless, reproduce the classical results obtained by virial expansion 
in the absence of short polymers. 

Various types of shape transformations of the long polymer can be obtained when the 
solvent acts as a good solvent for the long polymer and a poor solvent for short polymers and 
vice versa. These competitive interactions and the excluded-volume interaction cause the 
shape transformation of the long polymer. By systematically changing Flory’s y-parameter 
between the long polymer, solvent and short polymers, we investigate the shape transfor¬ 
mation of the long polymer induced by adding short polymers. 

This paper is organized as follows. In Sec. m the shape equation determining the optimal 
values of the expansion factor is derived from the partition function. We also present an 
alternative derivation of the shape equation from the mixing entropy using the effective 
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FIG. 1; Schematic illustration of the long polymer and short polymers, (a) Elongated coil state, 
and (b) compact globule state. R denotes the gyration radius of the long polymer, and dashed 
circles indicate the region occupied by the long polymer. Small red zigzags represent short polymers 
inside the region occupied by the long polymer. Small gray zigzags represent short polymers outside 
the region occupied by the long polymer. 

volume of the long polymer. In Sec. IIIIl we qualitatively examine the shape equation using 
an approximation. Numerical solutions of the shape equation are described in Sec. IIVI In 
Sec. E the chemical equilibrium between inside and outside of the region occupied by a 
long polymer is considered. Sec. IVll presents conclusion. 

II. DERIVATION OF THE SHAPE EQUATION 

In this section, we derive the shape equation of a long polymer from the partition function. 
We study the coil-globule transition of a long polymer with a degree of polymerization of nic 
in the presence of other polymers with a degree of polymerization of For both types of 
polymers, we introdnce the Knhn length of the monomer denoted by 1. [Ul We consider the 
case where rUc 3> mp. We focus on the coil-globnle transition of the long polymer and do 
not stndy the shape transition of short polymers here. The exclnded volnme of monomers 
of the long polymer is taken into acconnt in the partition fnnction. The short polymers are 
regarded as ideal chains. 
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A. Derivation by partition function 


The region occupied by a long polymer is characterized by the gyration radius R (Fig. 
d]). We also introduce the natural end-to-end distance of the long polymer denoted as Rq. 
Rn = t'y/2mc/3 can be obtained from the extremum of lnhFo(i?), where Wq{R) is given by 

HI 

W,(R) = fi^exp (-^) (1) 

for an ideal Gaussian chain. The expansion factor of the chain is given by a = R/Rq- The 
volume fraction occupied by the long polymer inside the region within R of the long polymer 
is expressed as, 

$ = m^iyR\ ( 2 ) 


is expressed using the expansion ratio as 

$ = 


ITlr 


(3) 


where s = (|)2 is a numerical factor independent of a and rric. 

The short polymers can move into the region occupied by the long polymer, as depicted 
in Fig. [H The number of short polymers inside the region occupied by a long polymer is 
denoted as n,. The volume fraction of short polymers inside the region occupied by the long 
polymer is denoted as mpPj and can be obtained from rUp'Pi = nirup^/R?. 


An expression for t 


le partition function in the absence of short polymers was described 


by Di Marzio. 19l-l21| In the presence of short polymers, the partition function can be 


generalized to (see Appendix for details) 

‘iR? 


W{R) = R^exp 


2mr P 


X 


Ui 




k=0 


'rric — l 

n 

. i=o 


l-ruppi-j— J 


X 


exp (1 - mpPi - $) - nicX^^pnipPi - mpniXp,s (1 - rUpPt - <!))] , (4) 


where and Xp,s represent Flory’s y—parameter between the long and short poly¬ 

mers, the long polymer and solvents, and the short polymer and solvents, respectively. 
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[10l | Flory’s y—parameters are given by the interaction energies divided by kBT, where 
ks is the Boltzmann constant and T is temperature. The o ccup ancy probability is as¬ 
sumed to be proportional to the volume fraction of emptiness. 19l-l2l| The factor given by 
(1 — mpPi — j£^/R^) represents the volume fraction of emptiness to place the j + lth segment 
of the long polymer when short polymers and j segments of the long polymer are already 
placed. Similarly, the factor given by (1 — krUpE'^/B?) represents the emptiness to place the 
k + 1th short polymer in the region occupied by the long polymer. 

In Eq. (jH), we rewrite the term relating the occupancy probability as 

me —1 / \ me —1 


^i.R) = n (^ ~ ^ 

3=0 




n 

3=0 




3 me 




(1 - iripPi) - j 


R 


w 


(1 - rripPi) + 1 


T (1 - rUpPi) -rric + lY 


(5) 


( 6 ) 


where r( 2 ;) is the gamma function. Using Stirling’s formula expressed by r(^-|- 1) {zjeY, 


22l | we obtain 
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f R\^ 

In n(i?) Ri ( j [(1 - rUpPi) In (1 - rUpPi) - (1 - rUpPi - d)) In (1 - rUpPi -$)]- rUc. (7) 


In the same way, we rewrite 

rii-l 


u{R) = ^ ~ 


k=0 




rii-l 


J-J- rUpi^ 


k=0 




R^ 


rUpP 


-k 


ITipP 


+ 1 




-rii + l 


( 8 ) 

( 9 ) 


Using Stirling’s formula, we obtain 

\nu3{R) 

From Eq. flTOl) . we hnd 

lna;(i?) 


Ui 


R^ 

rUpP 


(1 - rUpPi) In (1 - rUpPi) . 


( 10 ) 


p ^ —mpfii — 


R 


(1 - TUpPi) In (1 - rUpPi). 


( 11 ) 


We also applied Stirling’s formula to the rest of the factors in Eq. (jl]) to obtain 

3 


In 


R^yi 


riil \ P J 


piln{pi) + Ui. 


( 12 ) 
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The optimal shape of the long polymer can be obtained by maximizing \nW{R) with 
respect to a. We note that the volume fraction occupied by short polymers can be varied 
even if n* is not altered when R is changed. In this sense, rripPi depends on the volume 
fraction of the long polymer, Because Pi depends on a through $, it is convenient to 
define a new variable. 


r = 


rripPi mpTii 


(13) 


rric 

Here, rii can be regarded as a constant for the inhnitesimal small variation of R. In this 
work, we consider the small variation of R and regard r as being independent of a. 

By taking into account Eqs. l|7 |) - ffT^ . we obtain 


W{R) 2 
- TT- ~ a exp 


3/2 


—a — 


rric 3 


< Pi In {pi) + 1 - 


s(l + r) 


nira'^ 


In 1 - 


s(l + r) 




+ 




rric - rric X^,s + rXp.s + l+ r - 


TTlr 


where the effective y-parameter is defined by 

Xeff = - rAx + r'^Xp,s 


, (14) 


(15) 


and Ay is dehned by 


Ax X^,p Xp,s- 


(16) 


The optimal shape of the long polymer can be obtained by maximizing Eq. ffTTj) with 
respect to a. By calculating [a^/2]{d\nW{R))/{da) = 0 using Eq. fITT)) . we find 


a 


+ 


Srric 


1 + r 


r \ 3 

— + xXeffV^'5 + 

/ 2 


Ys 


In 


s(l + r) 


mcQc" 


= 0. (17) 


Equation ffTTl) was derived previously when r = 0 usIm the partition function except for a 


minor difference associated with the definition of a. 
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By further expanding Eq. flTTD in terms of s/{y/rn^a^), we obtain 


a' - 1 + 


3 rricr 


2 rrir. 


a^- 


3s 


rric [(1 + rf - 2yefr] “ ^ (1 + 


(18) 


For simplicity, we have expressed the volume of the polymer segment as £^. In fact, the 
volume of a polymer segment can be given by d‘^£ using the diameter of the polymer d. In 
this case, s is proportional to (R!£?. Below, we regard s as a parameter that characterizes 
the anisotropy of a polymer segment. 
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B. Derivation by free energy 


We further show that the partition function given by Eq. flT^ can be derived by combining 
the mixing entropy and entropic elastic energy of the long polymer when the mixing entropy 
is multiplied by the properly dehned volume occupied by the long polymer. 

Using the expansion factor a, the entropic elastic energy of the long polymer is given by, 

HI 

-^ = 0^-2 In a. (19) 

kbT 

We consider conhgurations of polymers using a lattice model. In the mean held approxi¬ 
mation, the mixing entropy at a site inside the long polymer is given by, jiu | 


5'mix (*h. Pi) = -kB \pi In (rUpPi) (1 - $ - nipPi) In (1 - d* - nipPi)] . 


( 20 ) 


In the absence of the elastic energy given by Eq, ([T9]), the free energy at a site inside the 
long polymer can be written as, 111 


-^in (^7 Pi) ^ 

+ X^,p'mpPi^ + x^,s (1 - mpPi - $) d* -t- Xp,s (1 - mpPi - $) rUpPi. 

Kb-L 

( 21 ) 

The entropic elastic energy for a long polymer is given in Eq. (|T^ and the free energy per 
site is given by Eq. fl2T]) . To add these terms, we should multiply the free energy of a site 
by a certain volume. The real polymer volume can be approximated by rricd'^i. We assumed 
that the effective volume occupied by a long polymer is 47tR^/{3£^). By comparing the site 
free energy given by Eq. fl2Tl) and the probability distribution given by Eq. flTTl) using the 
polymer volume given by AttR^/ [7| the total free energy Tin47ri?^/(3£^) including the 
elastic energy given by Eq. flT^ can be expressed as, 

^ ^ [Pi In {nipPi) -1- (1 - d) - rripPi) In (1 - d* - rripPi) -f x^,pmpPi^+ 


X^,s^ (1 - mpPi - <h) -t- Xp,s (1 - rripPi - mpPi] + 


kBT 


( 22 ) 


Equation is equal to — In TU (R) obtained from Eq. flTTD except for the unimportant 
term 1 -|- r — {r/nip). 


In the limit of pi = 0, Eq. (|22|) is simplihed to the equation obtained before. [2^ l27| 
The shape equation given by Eq. ffT7|l can be obtained by minimizing Eq. fl22|) with respect 
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to a. Previously, the equation representing the osmotic pressure balance was derived from 


Eq. fl2^ and the osmotic pressure balance was imposed by regarding the region occupied 
by the long polymer as being in a phase different from that of the other region, We did 
not impose the osmotic pressure balance at the boundary of the region occupied by the long 
polymer. Instead, we derive the optimal shape by minimization of Eq. fl22p with respect 
to the expansion factor. In the absence of short polymers, the resultant shape equation 


simplihes to the conventional form. 
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III. THE SHAPE EQUATION 


First, we study the shape equation derived from the partition function analytically. For 
convenience, we express Eq. flTSl) as 


/ (a) = 

where parameters are given by. 


, 3 

1 H— n, 
2 


3s 


a 


rricCtT 


^ = 0 




(23) 


Ct = ^, T = (l + r)2 
2/=y(l + r)^ 


2Xeff) 


(24) 

(25) 


and rii = rricr/mp from Eq. flT^ . When r = 0, Eq. fl23|) is simplihed to the same form of 
equation as that derived using the virial expansion. The term proportional to and 
in Eq. fl25]l represents the results of excluded-volume interactions, and that proportional to 
1/a^ represents the results of third-body interactions. 231-125 1 

Below, we summarize known results for Eq. fl23l) when Uj = 0 and generalize them for 


n* > 0. The results can be used when we study the coil-globule transition by solving Eq. 
flT7)) numerically. 


A. Extended state 

Extended states can be characterized by a > 1. If Ctr is positive, extended states are 
stabilized according to Eq. (1^ . When rric 3> 1, the extended states obey the scaling form 
~ ^/rn^CtT. By substituting s ~ and a ^ R/ into the scaling form, we 
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( 26 ) 


obtain the scaling relation for R in good solvents, 





when r>0. Ifr = 0,-R~ rric can be derived, which represents the same scaling form of 
the free polymers withont taking into acconnt the exclnded-volnme interactions. 


B. Compact State 


Compact states are characterized by a <C 1. When me 3> 1, the compact states obey 

m 


scaling form, — y'm^ctT ~ y/<y^- The scaling relation for R is obtained as H 


(27) 


when r < 0. The scaling of Eq. (l27ll represents the reasonable assnmption that the molecnlar 
weight of the long polymer is proportional to its volnme given by R? in the compact state. 


C. Coexistence state 

For the valnes other than a 3> 1 and a <C 1, we stndied Eq. (12^ analytically by assnming 
that CtT can take any valne and vary independently of y. The resnlts can be nsed as a gnide 
for nnmerical solntions of Eq. (na, where we take into acconnt proper constraints on CtX 
and y, in particnlar the fact that they are not independent fnnetions of the concentration of 
short polymers Pi. 

It is convenient to dehne 

g[a) = a^ - (^1 + (28) 

Eqnation fl2^ can be rewritten nsing g{a) as 

g{a) = y/m^CtT. (29) 

If g{a) has a maximnm and a minimnm in the region of a > 0, the coexistence of coil 
and globule states is possible depending on the value of r. The above condition can be 
restated as h{a) = g'{a)/a‘^ = 0 has two solutions for a > 0. The explicit condition can be 
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obtained by defining which satisfies both h'(a^) = 0 and h{a^) < 0 for > 0. Using 
Q/t = 2^ the condition is given by 0 < y < where 


yt = 


36 


44 X 53 


1 + = 


729 

32000 


, 3 

1 + jn, 


0.0228 I 1 + -ni 


(30) 


Therefore, the coexistence of coil and globnle states may be possible when 0 < y < 
depending on the right-hand side of Eq. ([291). 

We shonld also note that g{a) < 0 can be transformed into — 1 — (3/2)?7,j < y/a^. If 
0 < a < 1, he., — 1 < 0, g{a) < 0 is always satisfied becanse y is positive according to 
Eq. fl25|l . In addition, r < 0 is reqnired if the solntion of g{a) = ^/rn^CtT satisfies 0 < a < 1 
becanse g{a) = ^rricCiT < 0 and Ct > 0. The coexistence of coil and globnle states implies 
that we have two solntions and one of them should satisfy 0 < a < 1. Therefore, coexisting 
states can be obtained when r < 0. 


IV. SHAPE TRANSITION 

We next present numerical results of the shape transitions of the long polymer as a 
function of r. r is proportional to the number of short polymers inside the region occupied 
by the long polymer Uj. 

A. Parameter estimation 

We first estimate appropriate parameters to discuss the shape transition by considering 
r as a function of r (Fig. |2]). In general, the long polymer is in the elongated state if r > 0 
and in the compact state if r < 0, and the sign of r can be changed by changing r. By using 
Eq. flTHl) . Eq. fl2Tll can be rewritten as 

r = 1 - 2x$,, + 2 (1 + Ay) r + (1 - 2xp,s) (31) 

At r = 0, the elongated coil state can be obtained when < 1/2, and the compact globule 
state can be obtained when > 1/2. As for Xp,si when Xp,s > 1/2, t becomes negative 
upon increasing r, as can be understood from Eq. ([3T]). In contrast, when Xp,s < 1/2, 
r becomes positive upon increasing r. Based on the above reasoning, if * < 1/2 and 
Xp,s > 1/2 (Fig. |2] (a)), the coil-to-globule transition can be obtained, and if > 1/2 
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r 

FIG. 2: T as a function of r. r is proportional to the number of short polymers inside the region 
occupied by the long polymer, (a) Xp,s = 0.9, = 0.4, and x<i>,p = 0.18 (thick solid line), (b) 

Xp,s = 0.3, X'J>,s = 0.9, and x$,p = 0.433 (thin solid line), (c) Xp,s = 0.5, x$,s = 0.5, and x$,p = 0.1 
(dashed line). 

and Xp,s < 1/2 (Fig. [2](b)), the globule-to-coil transition occurs. Figure |2](c) presents the 
results when x$,s = 1/2 and Xp,s = 1/2. In this case, r slowly increases as r is increased. 
Below, we assume that rric = 1000, nip = 100, and s = 0.3387. 

B. Results 

1. Coil-to-globule transition 

Figure |3] shows the expansion factor a as a function of r with Xp,s > 1/2 and x<i>,s < 1/2 
where the coil-to-globule transition can be expected. In the absence of short polymers 
(r = 0), the long polymer extends in the good solvent characterized by x$,s < 1/2. The 
short polymers tend to aggregate in the poor solvent characterized by Xp,s > 1/2. When there 
are excluded-volume interactions between the long and short polymers, the long polymer 
may shrink if the short polymers tend to aggregate and the motion of short polymers is 
hindered by the presence of the long polymer. Therefore, the globule state can be obtained 
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r 


FIG. 3: a as a function of r when Xp,s = 0.9, X4>,s = 0.4, and X 4 >,p = 0.18. Here, r is proportional to 
the number of short polymers inside the region occupied by the long polymer. The red dashed lines 
were obtained using Eq. (11711 and the solid lines were obtained using the approximate expression 
given by Eq. (|M]l using nic = 10^, rup = 100, and s = 0.3387. 

when r exceeds a certain value. 

Figure |3] shows the coil-to-globule transition, a has two different values that correspond 
to the coil and globule states when r exceeds a certain value, and a single value (a > 1) 
when r is less than a certain value. Such a discontinuous transition suggests that it is a 
hrst-order transition and the compact globule state can be found when r exceeds a certain 
value. 

Incidentally, in Fig. [3l the dashed line was obtained using the original equation given by 
Eq. flT7|) and the solid line was obtained using the approximate expression given by Eq. fl23|l . 
We note that the difference between these lines is small and the qualitative r dependence is 
well captured by Eq. fl2^ . 

2. Globule-to-coil transition 

Figure H] shows the results for Xp,s < 1/2 and > 1/2 where the globule-to-coil 
transition can be expected. In the absence of short polymers (r = 0), the long polymer is 
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FIG. 4: a as a function of r when Xp,s = 0.3, X 4 >,s = 0.9, and x^,p = 0.433. r is proportional to 
the number of short polymers inside the region occupied by the long polymer. The red dashed lines 
were obtained using Eq. (11711 and the solid lines were obtained using the approximate expression 
given by Eq. (|M]l using nic = 10^, rup = 100, and s = 0.3387. 

in a compact state in a poor solvent characterized by > 1/2. In contrast, the short 
polymers are in good solvents. If a sufficient number of short polymers tend to be dispersed 
and involved in excluded-volume interactions with the long polymer, the extended coil state 
is expected. Figure 0] clearly shows that the long polymer changes from the compact globule 
state to the extended coil one as r increases. The transition is discontinuous. It should be 
noted that the extended coil state can be expected if a sufficient number of short polymers 
are involved in excluded-volume interactions with the long polymer regardless of the values 
of Xp,s and x^,p- Indeed, we obtain a similar hgure even when the short polymers can be 
regarded as a poor solvent for the long polymer characterized by = 0.6 > 1/2 if the 
other values in Fig. 0] are hxed. The discontinuous transition to the extended coil state 
is also obtained when Xp,s is changed to Xp,s = 0.6 and the other values in Fig. 0] are 
hxed. In the latter two cases, the transition to the extended state is mainly caused by the 
excluded-volume interactions between the long and short polymers. 
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FIG. 5: a as a function of r when Xp,s = 0.5, X 4 >,s = 0.5, and x^,p = 0.1. r is proportional to the 
number of short polymers inside the region occupied by the long polymer. The red dashed line was 
obtained using Eq. (fTTll and the solid line was obtained using the approximate expression given 
by Eq. (f23ll using rric = 10^, mp = 100, and s = 0.3387. 

3. Continuous coil expansion 

Figure |5] shows the results for Xp,s = 0.5, x<i>,s = 0.5 and = 0.1. The extended coil 
state expands gradually and the discontinuous transition is not obtained. 


Discussion 


In experiments, coil-to-globule transitions of DNA were observed when the so 


vent acted 


as a good solvent for both long DNA and short polyethylene glycol (PEG). 0, Although 
a coil-to-globule transition can be obtained according to our theoretical approach when 
the solvent acts as a good solvent for the long DNA, the condition for short polymers is 
different. By including the translational entropy of small ions, experimental results have 
been reproduced using mean held theory. Q Although such an ion effect is ignored in our 
theory, we systematically varied Flory’s y-parameters and obtained the general conditions for 
transitions between coil and globule states for neutral polymers. In some experiments, a coil- 


lb 






globule-coil transition was observed. j7|] Further study is needed to understand theoretically 
the globule-to-coil reentrance transition. 


V. EQUILIBRIUM BETWEEN OUTSIDE AND INSIDE THE AREA OCCUPIED 
BY A LONG POLYMER 


Here we discuss a way to determine the number of short polymers rii inside the region 
occupied by the long polymer. We consider the situation where short polymers can enter 
the region occupied by the long polymer and the concentration of short polymers inside the 
region occupied by the long polymer is different from that in the region outside the long 
polymer. The concentrations of short polymers outside and inside the region occupied by 
the long polymer can be expressed as and respectively. The relationship between 

the concentrations of short polymers in these regions can be determined by imposing the 
equality of chemical potentials. The free energy at a site outside the region occupied by the 
long polymer can be given by, [i^ 


Font 

k^T 


pin {rripp) + (1 — rUpP — rricc) In (1 — rUpP — rricC) + cln {rricc) + 


Xp,smpP (1 - TTipP - rricc) + x^,smcC (1 - rUpP - rrioc) + x#,pmcC nipP (32) 


in the mean held approximation, where c denotes the concentration of the long polymers 
outside the region occupied by a single long polymer. When short polymers can move across 
the boundary of the region occupied by a long polymer, the chemical potential outside the 
region containing the long polymer, p = dFout/dp, should be equal to that inside the area 
with the long polymer, dFin/dpi, which was determined using Eqs. (|2nil - (l2T|l . The result 
can be written as. 


In (p/pi 


rup In 


1 — rricC — TTipP 


1 - <F - rUpPi 

where Ay is dehned in Eq. (1161) . Equation 


= ^X^p ($ - moc) + 2xp,sml (p - pi ), (33) 


can be expressed as 


TTlpP 


(1 


nirC 


rrippY 


■ exp [AxmpmcC - 2yp,^mjp] = 


(1 _ ^ -mpPi)^^ [Aymp$ - 2yp,,mjp,] . (34) 

c and p for the regions outside and $ and pi for those inside the area occupied by a specihc 
long polymer are related by Eq. fl3T)) . 
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p 


FIG. 6: c as a function of p obtained using Eq. (f35]l when r = 0.5 and a = 2.90896, with Xp,s = 0.9, 
X$,s = 0.4, x^,p = 0.18, rric = 10^, rup = 100, and s = 0.3387. 


It is more convenient to rewrite Eq. flM|) in terms of r and a, where Pi is replaced by r 
using Eq. flT3|) . mpPi = r<I>. Equation flM|) can therefore be expressed as 



-AxmpUicC - {2xp,sr - Ax) rrip^^^ 
(l-(l + r)s/ 


mppexp {-2xp,smlp) 

(1 — nipP — rricc)"^^ 


(35) 


The optimal expansion factor, a, can be determined from Eq. (1T7|1 and Eq. (j35|l . 

In Fig. [HI we show c as a function of p for the parameters used to draw Fig. [3l The results 
indicate that the same values of a and r can be obtained if c and p satisfy the relationship 
represented by the line in Fig. [6l 


VI. CONCLUSION 

We theoretically investigated the excluded-volume and interaction effects of added poly¬ 
mer solvents on the transition between the coil and globule states of a long polymer. The 
excluded-volume effect was taken into account by generalizing the partition function of the 
long polymer to include the excluded-volume effect of short polymers. The interactions be¬ 
tween the segments of short polymers, those of long polymer and solvents were described 
using Flory’s y-parameters as in other theories. Following this approach, we obtained the 
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optimal expansion factors denoted by a from the maxima of the partition fnnction. To be 
more precise, the partition function was maximized by assuming that the number of short 
polymers n* inside the region occupied by the long polymer was conserved against the in- 
hnitesimal change in the gyration radius R of the long polymer. The shape equation was 
given in terms of the effective Flory’s y-parameter denoted by Xes- Xes was expressed by n* 
as well as Flory’s intrinsic y-parameters. 

The shape equation was simplihed using asymptotic expansion in terms of a. In the 
absence of short polymers, the simplihed equation reduces to the classical expression de¬ 
rived by the virial expansion. The simplihed shape equation captured the main features of 
the nonlinear shape equation and thus is very useful to analytically investigate the shape 
transition of polymers. Previously, the shape transition had been extensively studied using 
the simplihed equation obtained by the virial expansion in the absence of short polymers. 
Guided by this simplihed equation, we systematically studied the shape transition of a 


long polymer induced by the presence of short polymers. 

Using the simplihed shape equation and the ehective y-parameter XeS, we showed that 
coil-to-globule transitions can be induced by increasing n* when the solvent acts as a good 
solvent for the long polymer and a poor solvent for the short polymers. The long polymer 
tends to extend in a good solvent and short polymers tend to aggregate in a poor solvent. 
Under the competitive ehect induced by such a solvent, the long polymer shrinks if the long 
polymer segments are involved in excluded-volume interactions with a sufficient number of 
short polymers. These transitions were conhrmed by numerical calculations using the com¬ 
plete shape equation. In contrast, globule-to-coil transitions can be induced by increasing n* 
when the solvent acts as a poor solvent for the long polymer. This is because long polymer 
shrinks in the poor solvent, but expands through the excluded-volume interaction between 
the long polymer segments and short polymers. The globule-to-coil transition occurs at a 
lower concentration of short polymers when the solvent acts as a good solvent for the short 
polymers compared to that in the case where the solvent acts as a poor solvent for the short 
polymers. The globule-to-coil transition of a long polymer results from excluded-volume in¬ 
teraction with a sufficient number of short polymers and the threshold value for the number 
of short polymers can be low when the short polymers tend to be well dispersed in the good 
solvent. 

The excluded-volume effects can be also taken into account by considering the site-mixing 
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entropy of the lattice model. However, there was some ambiguity in dehning the effective 
volume of a long polymer to combine the site-mixing entropy and entropic elastic energy of a 
long polymer. We show that the minimization of the total free energy leads to the equation 
for the optimal expansion factor derived using the partition function if the mean volume is 
dehned by 47ri?^/(3£^), where £ is the Kuhn length and R denotes the radius of gyration of 
the long polymer. 



factor for unfolding states was expressed using the unoccupied volume to place the protein 
and the fraction of the reduced protein conformations due to prohibiting overlapping with 
crowding objects. By calculating the folding free energy using a reduction factor, the folding 
stability was discussed. Although the method of calculation is different, entropy reduction 
by excluded volume effect was taken into account both in the present theory and in the 
theory of protein folding stability. In this sense, this study gives another insight into the 
effect of macromolecular crowding using a single polymer chain. 

In this paper, we have implicitly assumed that the long polymer is flexible. Flexibility of 
conjugated polymers can be increased by formation of defects. Although conjugated poly¬ 
mers are stiff, various compact shapes of conjugated polymers were simulated by introducing 


defects. [37[ It has been pointed out that the folding characteristics of a long semiflexible 
chain differ from those of a short stiff chain. 38| Recently, it was demonstrated that a 
polymer of intermediate length formed various types of collapsed states, while the spherical 
collapsed state was obtained for l ong po lymers. 39| Spherical symmetric globules can be 


also classihed into several states. 


40|-|43| It was also shown that combined effect of rigidity 


and hydrophobicity is correlated with the folding kinetics and structures. 4^, 


451 To identify 


the shape of a long polymer, its stiffness, chemical interactions and hnite length should 
be explicitly considered. Although such complexity was ignored in this paper, the com¬ 
petitive and cooperative effects of excluded volume and the interaction between monomers 
and solvents were taken into account in this work using the partition function in the mean 
held approximation. By maximizing the partition function, we showed that various kinds of 
transition between the coil and globule states are possible by adding two different solvents. 


19 





















even if one of the solvents is short polymers. In particular, we qualitatively show that the 
shape transition can be induced by the excluded volume interactions if solvent acts as good 
(poor) solvent for a long polymer and poor (good) solvent for short macromolecules. The 
mechanism will apply for real polymers having stiffness and other chemical interactions. 
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Appendix A: Derivation of Eq. ([4]) 


If the excluded volume interactions are ignored, the total number of conhgurations of a 
long chain is proportional to i?^exp 2 m^)- presence of excluded volume interac¬ 

tions, the number of conhgurations is inhuenced by the presence of short polymers. 

We hrst consider the possible ways to place short polymers under excluded volume effect. 
When the hrst segment of the k + 1th short polymer is placed inside the region occupied by 
the long polymer, the number of possible ways to place the hrst segment is proportional to 
the free volume — krup, where krup denotes the total segments previously placed in 

the region. The second segment of the k + 1th short polymer can be placed to the nearest 
void if the space is allowed. When k short polymers have been previously placed randomly 
inside the area occupied by the long polymer, the probability of hnding a sufficient size of 
nearest void is given by 1 — kmp{£^/R^) in the mean held approximation. The probability of 
adding rup — 1 segments of the k + 1th short polymer after placing the hrst segment can be 
given by [1 — kmp{£'^/. In this way, the number of distinct conhgurations to place 
all Hi short polymers inside the region occupied by a long polymer is proportional to the 
factor 


'rii-l 

n 

.k=0 



^1 — knip 




(Al) 


where l/nj is multiplied to account for n^! ways of ordering indistinguishable short polymers. 

Now, we consider the excluded volume ehect of the long polymer segments under the 
presence of short polymers. For each conhguration whose incidence is proportional to 
R^ exp ( —requires that all rric successive segments of the long polymer can be 
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placed in a void. The probability of appearance of such a void can be estimated as follows. 
When short polymers and j segments of the long polymer are already placed, the next seg¬ 
ment of the long polymer can be placed in the nearest void and the probability of hnding 
the void is given by the volume fraction expressed by 1 — mpPi — in the mean held 

approximation. The probability of placing rric segments of the long polymer in a void can 
be given by [1 — nipPi — j{£^/R^)]. By taking into account all the above factors, we 

obtain 


np. 


rii-l 

n 

. k=0 


R^\ 


1 — km. 




'rric — l 

nil 

. j=0 


mpPi - j 


R^ 


(A2) 


By also taking into account the interaction energies and expansion potential of the long 
polymer, we obtain Eq. (|1]). 
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